Formation of a molecular Bose-Einstein condensate and an entangled atomic gas by 

Feshbach resonance 
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Processes of association in an atomic Bose-Einstein condensate, and dissociation of the resulting 
molecular condensate, due to Feshbach resonance in a time-dependent magnetic field, are analyzed 
incorporating non-mean-field quantum corrections and inelastic collisions. Calculations for the Na 
atomic condensate demonstrate that there exist optimal conditions under which about 80% of the 
atomic population can be converted to a relatively long-lived molecular condensate (with lifetimes 
of 10 ms and more). Entangled atoms in two-mode squeezed states (with noise reduction of about 
30 dB) may also be formed by molecular dissociation. A gas of atoms in squeezed or entangled 
states can have applications in quantum computing, communications, and measurements. 
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Introduction. — The recently discovered Bose- 
Einstein condensates (BEC), or matter waves, resemble 
in certain ways coherent electromagnetic radiation. This 
similarity stimulated the development of atomic optics 
0, involving non-classical states of the atomic fields, 
such as squeezed and entangled states Squeezed 
states are characterized by noise reduction, and can 
be applied in communications and measurements. 
Entangled states of a decomposable system cannot be 
expressed as a product of the component states, and can 
be used in quantum computing and communications. 
Squeezed atomic states can be formed in four-wave 
mixing in arrays of atomic traps Q, in multimode 
condensates || £j, in the decay of unstable BEC Q], in 
collisions of BEC wavepackets 0, and as the outcome 
of Bogolubov fluctuations subject to stimulated light 
scattering j|] or Belyaev dumping The squeezing 
can be measured experimentally by using homodyne 
detection, analogous to the one used in quantum optics 
(see Ref. Q). The key component of this method — a 
beam splitter — already exists (see Ref. |Io|). 

The present work suggests the dissociation of molecu- 
lar BEC as a source of atom fields in two-mode squeezed 
states that are entangled. Formation of single-mode 
squeezed states by the same mechanism has been dis- 
cussed 



Rcfs. [11 



|12| , pL 3| - Formation of entangled 
atomic pairs in the dissociation of individual diatomic 
molecules has been considered in Ref. [ 14 . Other mecha- 



nisms of formation of entangled gases have been discussed 
inRefs. 

The molecular BEC required as the source of the en- 
tangled gas is interesting in its own right, although it has 
not been realized yet. Formation of a molecular BEC by 
direct cooling of molecular gases is obstructed by the ro- 
tational degrees of freedom. An alternative method is 
the association of atomic BEC |ll], A process of 

photoassociation |l5| , realized experimentally Jl7| , is ob- 
structed by spontaneous emission Q. We consider here 
the association of atoms in BEC by Feshbach resonance 



||l6| in a time-dependent magnetic field. Such a process 
is associated with a large condensate loss observed in re- 
cent experiments Jll| . This loss follows from deactivation 
of the resonant molecules by inelastic collisions |ll| |o) , 
as well as from the formation of non-condensate atoms 
by molecular dissociation. An advantage of the use of 
Feshbach association is the possibility of reducing the 
negative effect of collisions by lowering the condensate 
density. Although the molecular BEC is made of excited 
molecules, and is therefore unstable, it still can be used 
as a source of atoms in entangled and squeezed states 
produced by the dissociation. 

Reference |2l[] treats the condensate loss as a dissoci- 
ation of single molecules. Many-body effects have been 
incorporated in Ref. [ p0[ by introducing a width to the 
molecular condensate state. A more rigorous analysis 
has been performed in Ref. p2| ], using second-order cor- 
relation functions. However, it was limited to the case 
of a time-independent zero detuning between the atomic 
and molecular states, and deactivating collisions were 
not taken into account. As in the present analysis, Ref. 
pof dealt with a time-dependent crossing of these states, 
in accordance with the experiments JS, and took into 
account deactivating collisions. We generalize here the 
parametric approximation, used in Refs. JlJ, for the 
description of molecular dissociation into a single atomic 
mode, and in Ref. |Q for a multimode analysis of fluctu- 
ations in an unstable BEC. 

The most impressive outcomes of the present approach 
are the extent of (near-total) conversion to entangled 
atoms or a molecular condensate, as well as the extreme 
degree of squeezing and the relatively long molecular 
BEC lifetimes achievable. 

The model. — Consider a system of coupled atomic 
and molecular fields (see Ref. ) described by annihila- 
tion operators in the momentum representation f „ (p, t) 
and ^ m (p,i), respectively. The coupling of the atomic 
and molecular fields (see Refs. ]l6|, ^(J) contains a prod- 
uct of two atomic creation operators and therefore de- 
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scribes the formation of entangled atomic pairs, in anal- 
ogy with parametric down-conversion in quantum optics 
(see Ref. ||). Spatial inhomogeneity due to the trapping 
potential and elastic collisions can be neglected here (see 
discussion below). Unlike the mean field used in Ref. 
p0{ , the atomic field is treated here as second-quantized, 
as in Ref. JT3| . 

Let the initial state of the atomic field at t = to be a 
coherent state of zero kinetic energy 



(p,i )|in) = (27r) 3/2 ^ ,5(p)|in 



(1) 



where \<po\ 2 = n a {to) is the initial atomic density and 
| in) is the time- independent state vector in the Heisen- 
berg representation. A pair of condensate atoms forms a 
molecule of zero kinetic energy. Therefore the resonant 
molecules can be represented by a mean field (p m (t) , such 
that 



(in|* m (p, t) in) = (2tt) 3/2 tp m (t) S (p) 



(2) 



where \ip m (t) \ 



(t) is the molecular condensate 



density. This approach (unlike Rcfs. (7[ |l3|) takes into ac- 
count the time dependence of the molecular mean field. 
Fluctuations of the molecular field due to collisions in- 
volving non-condensate atoms are neglected. 

The outcome of atom-molecule and molecule-molecule 
deactivating collisions is introduced, as in Ref. pOj , by 
adding molecular "dump" states. The elimination of 
these states in a second-quantized description is, how- 
ever, different. It is similar to the Heisenberg-Langevin 
formalism in quantum optics (see Ref. Q), but takes into 
account the nonlinearity of the collisional dumping. In 
the Markovian approximation, the equation of motion for 
the atomic field attains the form 

^* a (p, t) = H^ a (p, t) + 2g*<p m (t) (-p, t) + iF (p, f) 

(3) 

(using units with % = 1), where 



H 



P 

2m 



B (t) - S 



-ij\ip m (t) 



(4) 



and m is the atomic mass. The second term in H de- 
scribes the time-dependent Zeeman shift of the atom in 
an external magnetic field B (i), relative to half the en- 
ergy of the molecular state, which is chosen as the zero 
energy, fi is the difference in magnetic momenta of an 
atomic pair and a molecule, and Bq is the resonance value 
of B. The coupling of the atomic and the molecular fields 
g is related to the phenomenological resonance strength 
A as \g\ 2 = 2Tt\a a u\A/m (see Ref. [|(j)), where a a is the 
clastic scattering length. The parameter 7 describes the 
width of atomic states due to deactivating collisions (see 
Ref. [^0|). The quantum noise source F(p,t) provides 
conserving the correct commutation relations of the field 
operators. 



Generalizing the parametric approximation [il| , let 
us represent the atomic field operator in the form 



* a (p, t) = fi (p, t) i> c (p, t) + ^ (-p, t) ^ (p, t) 



x exp 



dti^\tp m (ti) 



(5) 



it, 



The operators i (p, t) are expressible in terms of 
ifr a (p, to), F (p, t), and the c-number solutions ip c>s (p, t) 
of the equations 

iip c>s (p, t) = Hi/j CtS (p, t) + 2g*(p m (t) ip* <c (p, t) , (6) 

given the initial conditions ip c (p, to) = 1, ip s (p, to) = 0. 
The atomic density 

n a (t) = (2vr)~ 3 J d 3 p 1 d 3 p 2 exp [i (p 2 - p x ) • r] 

x(in|*i(Pi.*)*a(P2,t)|m), (7) 

then appears to be r-independent, and comprises the 
sum n a (t) = no (t) + n s (t) of the densities of condensate 
atoms no (t) = (inl^a (0,t) |in)| 2 , and of non-condensate 
(entangled) atoms n s (t) in a wide spectrum of energies, 



n s (t) 



J dEn s (E, t) 



(8) 



where E is the kinetic energy of the non-condensate 
atoms. 

The equation of motion for the molecular mean field 
Lp m (t) is obtained by a similar elimination of the dump 
fields from the corresponding operator equation, followed 
by a mean-field averaging. We thus obtain 

iip m (t) = gm a (t) - i (jn a (t) + j m \(f m (t) | 2 ) ip m (t) , 

(9) 

where the parameter 7 m describes molecule-molecule de- 
activation collisions (see Ref. [^0|). Here m a is an anoma- 
lous density containing contributions of condensate and 
non-condensate atoms. The densities m a , no, and h s are 
all expressible in terms of ip c (p, t) and ij) s (p, t) . A numer- 
ical solution of Eqs. (||) on a grid of values of p, combined 
with Eq. (||), is consistently sufficient for elucidating the 
dynamics of the system. 

The present approach becomes mathematically equiv- 
alent to the approach of Ref. [^2| if the inelastic colli- 
sions are neglected and the detuning is time-independent. 
At a low molecular density, the effect of non-condensate 
atoms is equivalent to the contribution to the width of 
the molecular state made in Ref. [^0| for the same pro- 
cess. 

Formation of a molecular condensate. — Calculations 
were performed for two Feshbach resonances in colli- 
sions of Na atoms, using parameter values presented 
in Ref. (Eo|. The strong resonance, at 907 G, has the 



FIG. 1: Schematic illustration of transitions between atomic 
(AC) and molecular (M) condensates and non-condensed 
atoms (EA) on forward (FW) and backward (BW) sweeps. 
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FIG. 2: Time dependence of the densities of the atomic con- 
densate (dot-dashed line), the molecular condensate (dashed 
line), and the entangled atoms (solid line), calculated for the 
weak 853 G resonance in Na, with the initial atomic density 
no = 10 11 cm -3 , and the ramp speed B = —0.1 G/s (back- 
ward sweep). The dotted line shows the total atomic density 
(sum of the atomic densities and twice of the molecular one). 

strength A = 0.98 G, and the weak one, at 853 G, 
has the strength A = 9.5 mG. The difference of mag- 
netic momenta /i = 3.65 (in Bohr magnetons), the elas- 
tic scattering length a a = 3.4 nm, and the deactiva- 
tion parameters are 7 = 0.8 x 10~ 10 cm 3 /s and j m — 
10 _9 cm 3 /s. The neglection of elastic collisions is valid 
whenever no (to) <C 10 15 cm~ 3 for the weak resonance, 
and no (to) <C 10 17 cm -3 for the strong one. The spatial 
inhomogeneity can be neglected if the size of the conden- 
sate substantially exceeds 8 x 10 _2 cm -1 / 2 x n (to) 
and 2.5 x 10 _2 cm -1 / 2 x n (to), respectively, for the 
two resonances. Even when no (to) = 10 s cm -3 , these 
estimates set a minimal size of 8/im for the weak reso- 
nance and 2.5/mi for the strong one. The variation of the 
magnetic field is linear in time, B (t) = Bo + Bt. 

A relatively long-lived molecular condensate is formed 
more effectively in the case of a backward sweep, when 
the molecular state crosses the atomic one downwards 
(see Fig. [l]) , as proposed in Ref . [plj . The maximal con- 
version efficiency of the atomic condensate to a molecular 
one is 2max(n m ) /no ~ 0.8 for the weak resonance (see 
Fig. 0). On increase of the atomic density, or on dc- 
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FIG. 3: Conversion efficiency (solid line), optimal ramp speed 
B opt (dashed line), and lifetime of the molecular condensate 
Tm (dotted line) as a function of the initial atomic density, cal- 
culated for the weak resonance in Na (in a backward sweep). 

crease of the ramp speed, the conversion efficiency falls 
due to inelastic collisions. On increase of the ramp speed, 
the probability of crossing to the molecular condensate 
decreases, leaving more atoms in the atomic condensate 
(see Ref. p3[). At low atomic densities the conversion 
becomes less efficient due to a temporary gain of pop- 
ulation in the non-condensate atomic states. This ob- 
servation is peculiar to, and emphasizes the importance 
of, the simultaneous consideration of inelastic collisions 
and molecular dissociation within the second-quantized 
approach. (The approach used in J|(| does not describe 
the dissociation in the backward sweep.) 

Figure |j| shows that a substantial conversion efficiency 
is retained in a wide region of the condensate density, 
leaving much freedom in the choice of the ramp speed ap- 
propriate for experiments. The lesser the initial density, 
the longer the lifetime r m of the molecular condensate 
but the higher the precision required for the control of 
the magnetic field. 

The optimal ramp speed is approximately proportional 
to the initial density. This dependence minimizes the 
effect of a variation of parameters determining the con- 
version of the atomic condensate to the molecular one 
and loss of the molecular condensate. Indeed, the con- 
version to the molecular condensate is (in the fast de- 
cay approximation p0[) characterized by the parameter 
g 2 no/B. Similarly, the loss is characterized by the ra- 
tio of the deactivation lifetime (which is inversely pro- 
portional to the initial density), and the crossing time 
(which is inversely proportional to the ramp speed). 

The use of the strong resonance achieves a lower con- 
version efficiency, due to a gain in the temporary forma- 
tion of non-condensate atoms. The optimal ramp speed 
is more than two orders of magnitude larger than in the 
weak resonance, given the same initial density. 

Formation of an entangled gas. — As demonstrated 
in Ref. p|, the non-condensate atoms are formed in 
squeezed states, which now turn out to be two-mode 
squeezed states, as in Ref. Rj. It is similar to the state 
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two- mode squeezed states with the mean squeezing pa- 
rameter f fts 2.6, corresponding to a noise reduction of 
about 23 dB. The time dependence of the mean squeez- 
ing has a peak of f w 3.1 at t w —0.19 ms. The state of 
an entangled gas can be freezed at the peak time by fast 
turning off of the magnetic field. The energy spectra of 
the entangled-atom density and the squeezing parameter 
are presented in Fig. [|b. The density spectra are rather 
narrow, and the peak energy increases with time. The 
squeezing parameter reaches the value of r(E,t) w 3.5 
(noise reduction of 30 dB) at the energy E w 6 nK and 
the time t ~ 4 ms. 

Conclusions. — Both quantum corrections and deac- 
tivating collisions are necessary for the analysis of asso- 
ciation in atomic BEC due to Feshbach resonance in a 
time-dependent magnetic field. Over 80% of the atomic 
population can be converted to a molecular condensate in 
a backward sweep. The molecules dissociate onto atoms 
in two-mode squeezed states that are entangled. In a for- 
ward sweep, very high squeezing may be obtained with 
the parameter r reaching a value of more than 3. 
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FIG. 4: (a) Time dependence of the densities of the atomic 
condensate (dot-dashed line), entangled atoms (solid line), 
and the total atomic density (dotted line) calculated for 
the weak resonance in Na with the initial atomic density 
no = 10 14 cm -3 and ramp speed 50 G/s in a forward sweep. 
The dashed line shows the mean squeezing parameter f (t) 
[see Eq. (p"o|)]. (b) Energy spectra of the entangled-atom den- 
sity h a (E, t) (solid line) and the squeezing parameter r (E, t) 
(dashed line) calculated at the peak, t ~ —0.19 ms. The dot- 
dashed and dotted lines show their values on the plateau at 
t ^ 4 ms. 



of electromagnetic radiation formed by parametric down 
conversion. As in quantum optics Rcf. |^|, the amount 



of squeezing can he measured by the energy-dependent 



parameter r {E, t) related to the maximal and minimal 
uncertainties of quadratures involving creation and an- 
nihilation operators for opposite momenta Q]. A mean 
squeezing parameter, weighed by the spectral density of 
Eq. (|), 



F(t) = J dEh s (E,t)r(E, t) /r 



(*)■ 



(10) 



is used to describe the time variation of the squeezing. 

A stable gas of entangled atoms is formed by a forward 
sweep, in which the molecular state crosses the atomic 
one upwards (see Fig. |IJ). This process, too, is more ef- 
ficient in the weak resonance. The molecular density is 
then very low and persists a shorter time (compared to 
that in the backward sweep) due to fast dissociation. Fig- 
ure [|a demonstrates that more than 70% of the atomic 
condensate can be transformed into a gas of atoms in 
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